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Analyzing three approximate methods to locate liquid-solid coexistence in simple 



< 



a 



systems, an observation is made that all of them predict the same functional de- 
pendence of the temperature on density at freezing and melting of the conventional 



<£h Lennard-Jones system. The emerging equations can be written as T = Ap 4 + Bp 2 in 

o, 

^ ! normalized units. We suggest to determine the values of the coefficients A at freezing 



and melting from the high-temperature limit, governed by the inverse twelfth power 
repulsive potential. The coefficients B can be determined from the triple point pa- 



rameters of the LJ fluid. This produces freezing and melting equations which are 



exact in the high-temperature limit and at the triple point, and show remarkably 
good agreement with numerical simulation data in the intermediate region. 
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Several approximate approaches have been proposed to locate liquid-solid coexistence of 
different substances. This includes well known phenomenological criteria for freezing and 
melting, like e.g. Lindemann melting law, Hansen- Verlet freezing rule, Raveche-Mountain- 

X : 

$h ■ Street criterion for freezing, and a dynamical criterion for freezing in colloidal suspensions 

- - -' ft 

(for a review see Ref. []J). These criteria are typically based on the properties of only one 

of the two coexisting phases and predict quasi-universal values of certain structural or dy- 
namical quantities. Quasi-universality in this context means that a quantity is not exactly 
constant, but varies in a sufficiently narrow range for a broad variety of physical systems. 
Although these useful empirical rules can in most cases be easily implemented with a modest 
computational effort, sometimes it would be desirable to have simple analytical expressions 
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accurately describing phase coexistence, based only on the properties of the interparticle 
interaction potential. Even though such expressions can be non-universal (e.g. their appli- 
cability limited to a certain class of interactions), they can be quite helpful, especially in 
cases when interparticle interactions depend on a number of system parameters, which can 
vary from one situation to another «■ 

In this paper we analyze three such semi-empirical approaches developed to describe 
liquid-solid coexistence of simple systems and find that all of them predict the same func- 
tional dependence of the temperature on density at freezing and melting, when applied to 
the well known Lennard- Jones (LJ) potential of the form 

U(r) = 4e [(<x/r) 12 - [a/rf] , (1) 

where r is the separation between a pair of particles, and e and a are the energy and length 
scales. We therefore suggest to adopt this functional dependence and discuss how it can be 
used to produce accurate freezing and melting equations of the LJ system. The accuracy of 
the emerging equations is checked against existing results from numerical simulations. Some 
possible generalizations of this consideration are briefly discussed. 

As suggested in Ref. [3, for the potential U(r), which is a linear combination of repulsive 
potentials Ui(r), viz. U{r) = ^aj[/j(r), the freezing and melting curves are given, to a 
good approximation, by 

T L ' S (P) = $>7; L ' s (p), (2) 

where T i ' (p) is the temperature of the system of particles interacting via the potential 
Ui(r), as functions of the particle density p, along liquidus (L) and solidus (S). Equation (T2J) 
describes additivity of the melting curves.- Since the phase state of the system of particles in- 
teracting via the inverse-power-law (IPL) potential of the form U(r) = e(a/r) n is determined 
by a single dimensionless state variable (p<r 3 )(e/T) 3 / n , freezing and melting are specified by 
the value of the parameter c^' s = T L,s p _ri//3 , where T and p have been expressed in units of e 
and cr~ 3 , respectively- For the LJ system this immediately yields T L ' s (p) = 4c^ S p 4 — 4cg' s p 2 . 

Another approach is based on the cell model with a spherically averaged potential in 
the harmonic approximation. Assuming fee lattice in the solid phase and retaining only 
contributions from 12 nearest neighbors to the harmonic potential, the following melting 
equation has been derived in Ref. |5: 

T S = ^ 2 [*V(*)]', (3) 
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where x = (y/2/p) 1 ^ 3 and 5 is the Lindemann ratio. Applied to the LJ potential, Eq. (J3J) 
reduces to T s = 16<5 2 (llp 4 — 5p 2 ). The value of the parameter 5 can be chosen by fitting 
this equation to some known point on the melting curve. 

Finally, in a recent paper 6 it was observed that the freezing indicator in the form of 
the properly normalized second derivative of the interaction potential, C = [7"(A)A 2 /T 
remains practically constant along the freezing curve of the LJ fluid. Here A is the mean 
interparticle distance (related to the non-reduced density via A = p -1 / 3 ) and T is non- 
reduced temperature. The emerging equation for freezing of the LJ fluid can be written 
in reduced units as T L = (24/£)(26p 4 — 7p 2 ) The value of the constant C can be again 
determined by fitting this equation to some freezing point of the LJ fluid. For example, 
the high-temperature limit, governed by r -12 repulsion, can be used. In this way the value 
C ~ 290 has been estimated. 6 

All the considered approaches are approximate in their nature and involve serious simpli- 
fications . However, the very fact that three different (although not completely independent 
considerations) lead to the same functional form of the liquid-solid coexistence boundaries 
is indicative. It strongly suggests to use the functional form 

T L,S = ^L,S p 4 _ B L,S p 2 (4) 

to describe freezing and melting curves of the LJ system. Our main purpose here is not to 
discuss which of the approaches discussed above provides better values of the coefficients 
A and B. Instead, we suggest a procedure to determine the coefficients in Eq. (@J, which 
results in highly accurate freezing and melting equations of the LJ system applicable in the 
entire temperature range, from the triple point up to the high temperature limit. 

The coefficients A and B can be determined from the two reference points on the freez- 
ing/melting curve. One natural choice of the reference point is to take temperatures so 
high that the attractive contribution in the LJ potential can be ignored. In this regime 
the following dependence of density on temperature has been reported: p L ~ 0.813T 1//4 
at freezing and p s ~ 0.844T 1 / 4 at melting.-^ This fixes the values of the coefficients A at 
freezing and melting. We get A L o± 2.29 and A — 1.97. To determine the coefficients B 
in Eq. (J4]) it is appropriate to look at the low-temperature region of the liquid-solid coex- 
istence. In particular, triple point parameters of the LJ system are especially convenient 
for this purpose. The triple point temperature and densities obtained in different numerical 
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TABLE I: Triple point parameters of the LJ fluid obtained in different numerical simulations. 



T tT p\ p% Ref. 

0.660 0.862 0.960 [10] 

0.687 0.850 0.963 [11] 

0.692 0.847 0.962 [12] 

0.661 0.864 0.978 [13] 



simulations are summarized in Table [H There is some scattering of the data points, which is 
apparently due to differences in simulation methods as well as computational details such as 
system size, finite cutoff radius, etc. We will not discuss here the accuracy of each particular 
simulation. Instead, we take some "average" values for the triple point temperature and 
densities: T tr ~ 0.68, p\ t ~ 0.85 and p t s r ~ 0.96. These are also the values suggested by 
Hansen and McDonald. 9 This estimate is accurate to within few percent and can be further 
improved when consensus regarding the exact triple point parameters of the LJ system is 
reached. The emerging values of the coefficients B at freezing and melting are B L ~ 0.71 
and B s ~ 1.08, respectively. 

Expression (J4j) with the coefficients A and B derived above constitute freezing and melting 
equations of the LJ system, which are "exact" in the high-density high-temperature limit 
and at the triple point. In the intermediate regime, agreement with numerical simulation 
data is also remarkable, as shown in Fig. [T] The derived equations represent a considerable 
improvement over phenomenological fits proposed in the literature,— ~— which are sufficiently 
accurate only in the vicinity of the triple point. 

It is important to point out that similar freezing and melting equations can be constructed 
for a wide range of the LJ-type potentials. Obvious example is related to n-6 and more 
general Mie (n,m) model potentials. However, we do not elaborate further on this possibility, 
because triple point parameters of these systems are known with much lower accuracy.— 1 ^ 

To summarize, we propose very simple analytic equations for freezing and melting of the 
LJ system, which are very accurate in the entire range of temperatures and densities. These 
equations are superior to phenomenological fits proposed earlier. Similar equations can be 
easily derived for other model potentials of the LJ type. 
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FIG. 1: Comparison of the freezing and melting equations proposed in this paper (solid curves) 
with the representative set of numerical data related to liquid-solid coexistence of the Lennard- 
Jones system (symbols). For each value of the reduced temperature there are two values of reduced 
density. The left (right) symbol corresponds to freezing (melting). Stars, crosses, and triangles 



denote numerical results from Refs. 
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, and 
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